Zoosyst. Evol. 92 (2) 2016, 153-161 | DOI 10.3897/zse.92.8256 


> PENSUFT. 


museum fur naturkunde 


Life on the beach for a sand crab (Emerita rathbunae) (Decapoda, 
Hippidae): parasite-induced mortality of females in populations of 
the Pacific sand crab caused by Microphallus nicolli (Microphallidae) 


Juan Violante-Gonzalez', Scott Monks*, Guadalupe Quiterio-Rendon*, Sergio Garcia-Ibafiez’, 
Edvino Larumbe-Moran?, Agustin A. Rojas-Herrera! 


1 Unidad Académica de Ecologia Marina, Universidad Auténoma de Guerrero, Gran Via Tropical No. 20, Fracc. Las Playas, C.P. 39390, Acapulco, 


Guerrero, México 


2 Centro de Investigaciones Bioldgicas, Universidad Autonoma del Estado de Hidalgo, Apartado Postal I-10, C.P. 42001, Pachuca, Hidalgo, México 
3. Centro de Estudios Tecnologicos del Mar No 16, Interior Isla del Cayacal s/n, Zona Industrial, Cd. Lazaro Cardenas, C.P. 60950, Michoacan, México 


http://zoobank.org/499E 17E0-C26D-47D1-A479-02B5077849B6 


Corresponding author: Scott Monks (monks.scott@gmail.com) 


Received 24 February 2016 
Accepted 23 August 2016 
Published 8 September 2016 


Academic editor: 
Michael Ohl 


Key Words 


helminth parasite 

parasite-mediated mortality of hosts 
parasite population parameters 
Guerrero State 

Mexican coast 


Introduction 


Abstract 


Parasites, by definition, can affect mortality of their host, making parasitism an important 
biotic determinant of animal population dynamics and community structure. Reduction 
in the number of larger, reproductive age females in populations of the Pacific sand crab, 
Emerita rathbunae (Decapoda, Hippidae), was observed in studies of the helminth com- 
munity of this host. The aim of this study was to determine if high abundance of the meta- 
cercaria of the trematode, Microphallus nicolli (Microphallidae), causes mortality in this 
host. Females of £. rathbunae were collected from four sandy beaches in Guerrero State, 
Mexico, and helminths were collected from each crab. An analysis of variance (Anova) 
was applied to these data in order to identify differences in abundance between sizes of 
crabs, and an analysis of covariance (Ancova) was applied to identify differences in the 
abundance of metacercariae between locations. Parasite-related mortality was inferred by 
a decrease in abundance in older hosts. Linear and polynomial regressions of mean abun- 
dance of helminths (/og x+1 transformed data) vs. cephalothorax length of crabs were 
significant for the four populations of E. rathbunae, indicating increased mortality of 
older, more heavily infected female crabs and resultant removal from the population. En- 
capsulation and melanization of cysts by crabs was observed, indicating that an immune 
response by crabs also killed a portion of the cysts from subsequent exposures. Mortality 
of hosts through behavioral modification favoring transmission of highly infected crabs 
was suggested as the driving force behind this process. 


predation (Hansen and Poulin 2005; Koehler and Poulin 
2010; Latham and Poulin 2002c; Poulin 1999; Thomas 


Parasitism is an important biotic determinant of animal 
population dynamics and community structure. Parasites 
can influence host numbers within a population by either 
decreasing recruitment rate (by reducing host fecundity) 
or increasing mortality rate. Host mortality attributable 
to parasites can occur as a reduction in host survival due 
to the pathological consequences of parasitic infection, 
or as alterations in a host’s phenotype and/or behavior 
to facilitate parasite transmission to definitive hosts via 


et al. 1995). 

Several approaches based on statistical evidence have 
been used in a number of studies in an attempt to deter- 
mine if infection with parasites affects host mortality rates 
in natural populations (Koehler and Poulin 2010; Latham 
and Poulin 2002c; Thomas et al. 1995). A frequent fea- 
ture of many models is that host mortality rises as par- 
asite load increases (Koehler and Poulin 2010; Latham 
and Poulin 2002c). In a system where the host continues 
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Figure 1. Helminth parasites of Emerita rathbunae. A Microphallus nicolli, B Everted proboscis of an excysted specimen of Pro- 
filicollis sp.; C Cestode larva; D Proleptus sp., inset (*) = anterior end of worm. Scale bars: A = 500 um; B = 200 um; C = 100 um; 


D=100 um. 


to accumulate parasites over time, and assuming no den- 
sity-dependent regulation of parasite numbers, older indi- 
viduals can be expected to have higher average parasite 
loads than younger hosts. In some cases, however, older 
hosts have fewer parasites on average than intermediate 
age hosts, suggesting that more heavily infected individ- 
uals are removed from the host population by some effect 
of the parasites (Koehler and Poulin 2010; Latham and 
Poulin 2002c; Thomas et al. 1995). 

The Pacific sand crab, or mole crab, Emerita rathbu- 
nae Schmitt, 1935 is common in the splash zone of san- 
dy beaches in the intertidal zone (Hendrickx and Harvey 
1999; Rios-Eldsegui and Hendrickx 2015). It is a good 
model for parasitological studies due to its wide geographi- 
cal distribution, from the southeastern Gulf of Califor- 
nia, Mexico, to the southern coast of Ecuador (Tam et al. 
1996), and its importance in the macrofauna food chain of 
sandy beaches, where it can be very abundant (Oliva et al. 
2008). A number of studies on the parasites of E. analoga 
have been made (Alvitres et al. 1999; Iannacone et al. 
2007; Oliva et al. 2008; Oliva et al. 1992: Smith 2007), 
but only one study of the parasites of E. rathbunae (re- 
ported as E. analoga) has been made in Mexico (Violan- 
te-Gonzalez et al. 2012; Violante-Gonzalez et al. 2015). 
In that study, the authors (Violante-Gonzalez et al. 2012; 
Violante-Gonzalez et al. 2015) found that specimens 
of EF. rathbunae inhabiting the coasts of Guerrero were 
parasitized by four species of larval helminth: one me- 
tacercariae [Microphallus nicolli (Cable and Hunninen, 
1938)]; one cystacanth (Profilicollis sp.); one cestode 
(Trypanorhyncha); and one nematode (Proleptus sp.) 
(Fig. 1A—D). Infection levels of M. nicolli were signifi- 
cantly higher than reported for populations of FE. analoga 
(hundreds to thousands per infected crab) from another 
North American location (Smith 2007). This suggested 
that these helminths have the potential to exert a substan- 
tial effect on populations of E. rathbunae. The objective 
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of the present study was to test that hypothesis and de- 
termine if high abundance of M. nicolli in populations of 
E. rathbunae on the coasts of Guerrero, Mexico, effects 
mortality rates of this species of host in a manner reported 
for other species of crustacean at other locations. 


Materials and methods 


Sampling 


The Mexican Pacific coastal region experiences two distinct 
climatic seasons, a rainy period from June to November 
(precipitation ~ 430 mm), and a dry season from Decem- 
ber to May (precipitation < 70 mm). Adult females of E. 
rathbunae (Table 1) were collected from four sandy beaches 
on the Pacific Coast of Guerrero, Mexico, between August 
and December 2009: San Andrés (16°42’ N, 99°40’ W, n 
= 146); Revolcadero (16°47’ N, 99°47’ W, n = 114); Las 
Trancas (16°59’ N, 100°13’ W, n= 173); and Ixtapa (17°39’ 
N, 101°36’ W, n = 68) (Fig. 2). Males of E. rathbunae are 
much smaller than females, and measure from 2-14 mm 
(i.e., dwarf males), and they generally remain attached to the 
coxa of the female’s pleopods, or they are located among the 
eggs (Rios-Elosegui and Hendrickx 2015). Therefore, due 
to their small size they were not considered for this study. 
Female crabs were captured by hand, placed in labeled 
plastic bags and transported to the laboratory of the Uni- 
dad Académica de Ecologia Marina, Universidad Auténo- 
ma de Guerrero, Acapulco. Before dissection, cephalo- 
thorax length (CL in mm) was measured to the nearest 
0.1 millimeter using a digital caliper. This parameter was 
treated as the standard crab size in the analyses. Crabs 
were dissected by first removing the carapace, and then 
examining the internal structures, hepatopancreas, stom- 
ach, intestine, and muscle tissue, for helminths. Helminths 
were collected using standard practices and processed 


Zoosyst. Evol. 92 (2) 2016, 153-161 


155 


Table 1. Infection parameters of metacercariae of Microphallus nicolli in Emerita rathbunae from four coastal locations in Guer- 
rero, Mexico. Significant values for the relationship between crab size (cephalothorax length) and metacercariae abundance per 
infected host (r,) are in bold. P (%) = Prevalence of infection (% of infected crabs); Range of Intensity = minimum and maximum 


number of metacercariae; r, = Spearman’s correlation coefficient. 


: Cephalothorax | Total number Mean abundance of |_ Range of 
Date (2009) Locality No. of crabs 1 P (%) ' ; r 
length (mm) of parasites metacercariae Intensity : 
Sep. San Andrés 146 SiG ee.O 40,3414 100 2763.1 + 1235.6 | 300-6000 | 0.208 
Oct. Revolcadero 114 34.24 1.7 LO7§921 100 946.7 + 796.8 103 - 4235 | 0.609 
Nov. Las Trancas 173 Sole eee 403,177 100 2330.5+ 1076.8 | 1050-5800} 0.444 
Dec. Ixtapa 68 SAN E56 58,241 100 856.5 £59775 15-2550 | 0.689 


Guerrero 


Pacific 
Ocean 


Acapulco 


Figure 2. Map indicating the four collection localities along the 
Pacific coast of Guerrero. 1. San Andres; 2. Revolcadero; 3. Las 
Trancas; 4. Ixtapa. 


for identification (Pritchard and Kruse 1982). Infection 
parameters included prevalence (percent infected hosts); 
abundance (number of parasites per examined host) ex- 
pressed as the mean + standard deviation; range of intensi- 
ty (Bush et al. 1997); and volumetric abundance (volume 
of live metacercariae of all development stages per host, 
in mm>*) (see Violante-Gonzalez et al. 2015 for discussion 
of these data). Vouchers of specimens were deposited in 
the Coleccion de Helmintos, Universidad Autonoma del 
Estado de Hidalgo, Mexico (CHE-F-0019). 


Statistical analyses 


Infection parameters included prevalence (percent of in- 
fected hosts); abundance (number of parasites per exam- 


ined host) expressed as the mean + standard deviation; 
range of intensity (Bush et al. 1997); and volumetric 
abundance (volume of live metacercariae of all develop- 
mental stages per host, in mm*). An analysis of variance 
(Anova) was applied to identify differences in sizes of 
the crabs, and an analysis of covariance (Ancova) was 
applied to identify differences in the abundance of meta- 
cercariae between locations. Data were log-transformed 
when significant deviations from normality were iden- 
tified. Correlations were calculated using the Spearman 
range coefficient (r,). Significance for all the statistical 
analyses was established at P = 0.05, unless stated oth- 
erwise. The software SPSS (Version 20) was used for all 
statistical analyses. 

Size varied among the encysted metacercariae of M. 
nicolli, and a number of melanized cysts were found in 
each infected crab, so samples of cysts were measured 
(length and maximum width in wm) using an ocular mi- 
crometer for two locations (San Andrés and Las Trancas; 
n= 866). These measurements were used to calculate the 
volume of live cysts (mm?) using the volume formula for 
a scalene ellipsoid [V = (4/3) z (a/2) (6/2) (c/2) where 
a 1s cyst length, b is cyst width, and c is cyst height]. 
Melanized cysts were considered to contain dead meta- 
cercaria (Bryan-Walker et al. 2007; Koehler and Poulin 
2010). Cysts also were classified into four categories 
based on size and characteristics (live or dead/melanized), 
following a modified version of the classification report- 
ed by Keeney et al. (2007): 1. = immature cyst (small); 2. 
= late immature (medium); 3. = mature (large); and 4. = 
melanized (dead). Crabs from San Andrés and Las Tran- 
cas were pooled and divided arbitrarily into 6 equal (2.1 
mm) size-classes 1n order to describe cyst distribution by 
size (cephalothorax length = CL) of crabs. 

Mean parasite load (total number of live metacercari- 
ae of all sizes/total number of crabs that were examined) 
was plotted versus size-class of crabs for the four local- 
ities in order to evaluate the affect of this parameter on 
host mortality; a close relationship was assumed between 
age and size (Latham and Poulin 2002b). The variance to 
mean abundance ratio, a measure of parasite aggregation 
(Bush et al. 1997), was calculated to determine variation 
in infection levels by host size-class. Differences in mean 
parasite abundance and the variance to mean abundance 
ratio were analyzed as a function of increasing size class. 
Linear and polynomial (curvilinear) regressions were 
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then fitted using the log-transformed data (log x+1) of 
mean parasite abundance (dependent variable) and CL 
(independent variable), by location. The regressions pro- 
viding the best fit for the data were selected using the 
highest coefficient of determination (r?) value and a sub- 
stantially low P value for each regression (Latham and 
Poulin 2002b). 


Results 


The total sample size was 591 female crabs (E. rathbu- 
nae), ranging from 68 individuals at Ixtapa to 173 at Las 
Trancas (Table 1). Cephalothorax length was signifi- 
cantly different between collection locations: 32.9 mm 
+ 5.6 mm (Ixtapa) to 38.2 mm + 2.2 mm (Las Trancas) 
(Anova, F = 80.8, P < 0.0001) (Table 1). Encysted me- 
tacercariae (Fig. 3) of M. nicolli were recovered largely 
from the hepatopancreas of infected crabs. Prevalence of 
M. nicolli was 100% at all four locations, although mean 
abundance varied from 857 metacercariae/crab + 598 me- 
tacercariae (Ixtapa) to 2,763 metacercariae/crab + 1,236 
metacercariae (San Andrés, Table 1) (Ancova, F = 43.62, 
P=0.0001). 

Cephalothorax length was correlated positively with 
abundance of metacercariae at all locations (P < 0.05), 
although correlation coefficient values (r,) were low in 
most cases (Table 1). The overall percentage of live (un- 
melanized) cyst decreased as the developmental stage 
progressed from immature to mature cysts (Table 2). The 
percentage of melanized cysts corresponded to the three 
previous development stages, although the general cyst 


Figure 3. Cyst of an older metacercaria of Microphallus nicolli. 
Scale bar: 100 um 
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pattern remained unchanged, 1.e., all crabs had a higher 
percentage of immature than mature but melanized cysts. 

Distribution of cysts, based on developmental stage and 
crab size, did not differ at the two locations where this was 
analyzed (San Andrés and Las Trancas). The percentage 
of mature cysts decreased as crab size increased, and the 
percentage of melanized cysts increased with the increase 
(small to large) in the size of crabs (Fig. 4A, B). Mean 
volumetric abundance of metacercariae, considering the 
three developmental stages, decreased in the largest size 
classes of crabs at both locations (Fig. 5A, B). 

Linear and polynomial regressions of the relationship 
between mean abundance of helminths (log x+1) and 
cephalothorax length of crabs were significant at all four 
beaches (Table 3). However, the coefficients of determi- 
nation (1?) values for the polynomial regressions were 
higher, indicating that the curvilinear regressions provid- 
ed a better fit to the data. Mean abundance of M. nicolli 
plotted against crab size-class produced a convex curve 
in the San Andrés and Ixtapa populations, indicating that 
cyst number tended to decline in the largest size-classes 
(Fig. 6A, D). In contrast, the curves generated for Revol- 
cadero and Las Trancas were concave (Fig. 6B, C). The 
degree of parasite aggregation was lower in the largest 
size-classes at San Andrés and Ixtapa (Fig. 6A, D), which 
coincides with the convex curves produced by the poly- 
nomial regressions for these two locations. 


Discussion 


Species of Microphallus infect several species of marine 
crustaceans that function as second intermediate hosts; 
adult worms mature mainly in coastal birds (Hansen and 
Poulin 2005; Pina et al. 2011). The metacercariae exhib- 
it high infection site-specificity and encyst mainly in the 
hepatopancreas of crabs, where they reach high prev- 
alences (Hansen and Poulin 2005; Koehler and Poulin 
2010; Pina et al. 2011; Robaldo et al. 1999). In crusta- 
ceans, the hepatopancreas is the main organ involved in 
carbohydrate metabolism (glucose and glycogen), which 
means that it provides greater availability of energy re- 
sources (nutrients) for establishment of helminths (Ro- 
baldo et al. 1999). All (100%) of the specimens of E. 
rathbunae from the four locations were infected with M. 
nicolli, recovered mainly from the hepatopancreas, con- 
firming the expectation of a strong habitat specificity for 
members of this genus (Koehler et al. 2011). 

Although prevalence of M. nicolli did not differ be- 
tween the four locations, mean abundance was signifi- 
cantly higher at San Andrés and Las Trancas (Table 1) 
than at the other two locations. Temporal and spatial 
variations in levels of infection in intermediate hosts are 
commonly observed patterns in studies of parasite ecol- 
ogy (Latham and Poulin 2002b; Oliva et al. 2008; Smith 
2007; Studer and Poulin 2012). Spatial variation in the 
level of infection such as that observed here has been 
attributed by other authors to differences in abundances 
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Figures 4. Percentages of cysts of Microphallus nicolli by 
development stage of cysts and size-classes of female crabs 
(Emerita rathbunae). Hatched bars = immature cysts; Shaded 
bars = late immature cysts; Open bars = mature cysts; Dark bars 
= melanized cysts. Data are only from San Andrés (A) and Las 
Trancas beaches (B). Stages of cyst maturity were modified 
from Keeney et al. (2007). 
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Figure 5. Relationship between mean abundance of cysts of Mi- 
crophallus nicolli calculated from volumetric measurements and 
the cephalothorax length of females of Emerita rathbunae. Data 
are from only San Andrés (A) and Las Trancas beaches (B). 
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Figure 6. Relationship between mean abundance of cysts (log x + 1) of Microphallus nicolli and size-classes of females of Emerita rath- 
bunae based on cephalothorax length. Data are for all four localities: San Andrés (A), Revolcadero (B), Las Trancas (C), and Ixtapa (D). 
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Table 2. Mean length and width measurements and the volume (mm*) of encysted metacercariae of Microphallus nicolli by devel- 
opment stage. Global percentage of melanized cysts corresponds to the sum of the percentages found in the three stages. 


: Volume | Global percentage of | Mean number of 
Cyst stage No. of cysts Length (mm) Width (mm) ; : 
(mm?) melanized cysts melanized cysts 
Immature 247 0.259 + 0.040 0.221 + 0.034 0.0066 37% 91.4 
Late immature Za 0.350 + 0.038 Os Ole tO1035 0.0166 26% 86.2 
Mature L651 0.591 + 0.046 O1393'1;03038 0.0477 16% Bae 
Melanized cyst 0.317+0.062 | 0.269+0.044 | 0.0125 21% 


of final and intermediate host (lannacone et al. 2007; 
Latham and Poulin 2002b; Latham and Poulin 2002c; Ol- 
iva et al. 2008; Smith 2007). In this study area, a possible 
first intermediate host for M. nicolli is the snail Agaronia 
testacea (Lamarck, 1811), and the probable final hosts 
are coastal marine birds; species of Scolopacidae and 
Recurvirostridae (Charadriiformes); members of 7ringa 
and Numenius are known to actively feed on EF. rathbu- 
nae (Violante-Gonzalez et al. 2012; Violante-Gonzalez et 
al. 2015). However, to date, no studies have been made 
to identify precisely the definitive host of M. nicolli in 
Guerrero. 

It is possible that variation in the abundance of M 
nicolli between the four locations is the result of a fac- 
tor or process related to the size of the crabs that might 
influence this infection parameter (Poulin 1999). A pos- 
itive correlation between crab size and abundance was 
observed at all four locations (Table 1), with mean abun- 
dance highest in crabs from San Andrés and Las Tran- 
cas (Table 1). However, this tendency also is consistent 
with the hypothesis that largest crabs (1.e., older) have 
accumulated parasites for a longer period of time. Indi- 
viduals of E. analoga can to live up to 3.5 years (Oliva et 
al. 2008), so individuals of that species accumulate large 
quantities of helminths during their lifetime. In our study, 
it was not possible to determine the strength of a sim- 
ple age-infrapopulation size effect because the maximum 
age/size of E. rathbunae has not been documented, but 
our results are similar to the mortality seen in studies of 
E. analoga (Alvitres et al. 1999; Iannacone et al. 2007; 
Oliva et al. 2008; Oliva et al. 1992; Smith 2007) and oth- 
er shore crabs (Latham and Poulin 2002a; Latham and 
Poulin 2002b). 

Metacercaria of M. nicolli continue to grow post-en- 
cystment (Peoples and Poulin 2011; Pina et al. 2011), 
and we found that size of cysts differed according to 
developmental stage (Table 2). Mean volume of imma- 
ture (small) cysts was 0.0066 mm+°, and the volume of 
mature cysts was 0.0477 mm?, representing an increase 
from initial volume of up to 700% during encystment in- 
side the crabs. However, the mean number of cysts/crab 
decreased with developmental stage (247, 233, and 151, 
respectively). We interpret this as support for increased 
mortality of crabs as metacercariae develop. There are no 
studies of the histopathology caused by these helminths, 
but the increase in volume of metacercariae undoubtedly 
damages the tissues of the crabs. It also 1s not known if 
the size of metacercaria is linked directly to the size of 
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crabs (Ruiz-Daniels et al. 2013). These are areas for fur- 
ther research. 

A high percentage of cysts had been encapsulated 
and melanized (Table 2). This process, encapsulation 
and subsequent cellular melanization of metacercariae, 
is a common response of the host’s immune system in 
crustaceans, and it results in the death of the metacercaria 
(Keeney et al. 2007; Koehler and Poulin 2010). Immature 
cysts frequently were melanized more than were those of 
older stages, and the percentage of melanized immature 
cysts increased progressively from the smallest to the 
largest crab size classes (Fig. 4A, 4B). This suggests that 
the immunologic response in E. rathbunae is age-depen- 
dent; 1.e., through this response larger (older) crabs can 
eliminate a large percentage of recently-introduced para- 
sites, thus slowing greater accumulation of parasites and 
concurrent mortality (Keeney et al. 2007; Koehler and 
Poulin 2010). As well, the mean number of melanized 
cysts decreased with developmental stage; we cannot ex- 
plain this decrease in number and there is no available ev- 
idence that the melanized cysts eventually were absorbed 
by the crabs. There also is no information on the potential 
for genetically mixed populations of M. nicolli (brought 
to these beaches by birds from completely different re- 
gions), although this is likely. This could confound our 
results if there are different immune responses to different 
strains, as was seen by Koehler et al. (2011). Molecu- 
lar studies of individual helminths (DNA fingerprinting) 
should be carried out in the future to resolve this issue. 

The overall percentage of mature (large) cysts was 
lower than the percentage of the two earlier stages (Ta- 
ble 2), which could be attributed to larger cysts requiring 
more nutrients for growth, thus causing death of imma- 
ture stages, although there is no direct evidence for this. 
Intraspecific competition for resources at high cyst con- 
centrations (800 to 2000 cysts in the present study) 1n the 
same structure (hepatopancreas) could be causing cyst 
death. Density-dependent effects acting on regulation of 
metacercarie growth are common in marine crustacean 
hosts (Fredensborg and Poulin 2005; Keeney et al. 2007), 
although more study is necessary to elucidate direct and 
indirect causes. This regulator effect on survival of par- 
asites according to size may be very important for M. 
nicolli because, unlike many other microphallid larvae, it 
grows several hundred times its initial volume inside the 
second intermediate host (Fredensborg et al. 2004; Kee- 
ney et al. 2007; Koehler and Poulin 2010; Latham and 
Poulin 2002b). 
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Ample evidence exists of parasite-induced mortal- 
ity in marine crustacean populations under field condi- 
tions (Fredensborg et al. 2004; Koehler and Poulin 2010; 
Latham and Poulin 2002b; Thomas et al. 1995). Para- 
site-induced host mortality can be deduced by means of 
a decrease in infection intensity in older hosts, which we 
observe in our data. Using an acanthocephalan/crab sys- 
tem, Latham and Poulin (2002) suggested that reductions 
in mean parasite numbers in the largest crab size-classes 
could be attributed to the loss of heavily-infected individ- 
uals; this also was observed in our data. With a similar 
approach, Koehler and Poulin (2010) reported host mor- 
tality in four species of crab from New Zealand, and they 
concluded that this was induced by two microphallids 
(Maritrema novaezealandensis and Microphallus sp.). 

If an infrapopulation of a species of parasite does not 
induce host mortality, then a positive linear relationship 
can be expected between host size and parasite abundance 
due to parasite accumulation in the host over time (Koehler 
and Poulin 2010). In the data presented herein, polynomi- 
al (curvilinear) regressions had a better fit to mean abun- 
dance data per crab size-classes at all four locations (Table 
3) than linear regressions. However, the curves generated 
were of two types; concave and convex. For Revolcadero 
and Las Trancas, the best-fit curves were concave (Fig. 
6A, D), suggesting an increase in parasite load with crab 
size surpassing a linear relationship. The regression curves 
for San Andrés and Ixtapa were convex (Fig. 6B, C), in- 
dicating a reduction in parasite load in the largest (oldest) 
size Classes. The latter could be attributed directly to death 
of those hosts harboring greater numbers of parasites and 
removal of both host and parasites from the population of 
crabs. It is interesting that these differences do not appear 
to be related to parasite intensity (Table 1). 

In some studies of host mortality induced by helminth 
parasites, differences have been reported in the effect of 
parasites between host species and locations (Koehler 
and Poulin 2010; Latham and Poulin 2002b). These dif- 
ferences have been attributed to possible host manipula- 
tion by the parasites in order to increase the likelihood 
of predation by definitive hosts, and the presence of the 
final host (1.e., shorebirds) in locations where host popu- 
lation regulation was studied. Although several species 
of Microphallus can manipulate host behavior (Hansen 
and Poulin 2005; Ruiz-Daniels et al. 2013; Thomas et al. 
1995), there have been no studies of the possibility of in- 
fected individuals of E. rathbunae being manipulated by 
M. nicolli. There is, however, a high abundance of marine 
coastal birds at San Andres (unpublished data) that could 
cause heavy predation of larger, highly-infected crabs, 
and, during this study, migratory birds were observed 
catching and eating crabs on these beaches. This might 
explain the possible host population regulation suggested 
by the regression analysis (Table 3) and the mean para- 
site abundance vs. crab size class graph (Fig. 6A). This is 
not the case for Ixtapa, although the statistical analyses 
also indicated a reduction in the parasite load in the larg- 
est size classes of crab (Fig. 6D). The number of birds 
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Table 3. Results of the regression analysis used to evaluate the 
relationship between number of cysts of Microphallus nicolli 
(log x+1) per crab and the cephalothorax length of crabs from 
four coastal locations in Guerrero, Mexico. n = number of crabs 
that were examined; R? = coefficient of determination; P = sig- 
nificance level. Note: 100% of the crabs examined were infected. 


Linear Polynomial 
Locality n R? P R? P  |Curve type 
San Andrés | 146 | 0.297 | 0.0001 | 0.330) 0.0001] Convex 
Revolcadero | 114 | 0.393 | 0.0001 | 0.395 | 0.0001 |} Concave 
Las Trancas | 173 | 0.235 | 0.0001 | 0.300 | 0.0001 |} Concave 
Ixtapa 68 | 0.508 | 0.0001 | 0.588 | 0.0001] Convex 


visiting that beach noticeably was less, probably due to 
constant tourist traffic, which would directly affect the 
parameters of the helminth population in E. rathbunae. 

The most plausible explanation for the lower abun- 
dance of metacercaria in larger size crabs is the regula- 
tion of the population of E. rathbunae through elimina- 
tion of a high proportion of metacercariae cysts, directly 
by mortality of crabs or by favoring the capture of the 
more highly-infected crabs (the larger, older hosts) (Kee- 
ney et al. 2007; Koehler and Poulin 2010). As we ob- 
served, E. rathbunae can eliminate/inactivate a consid- 
erable percentage of cysts through melanization, most 
of which are immature. The largest crabs may have a 
stronger immune response, allowing them to eliminate 
a greater number of cysts by encapsulation, followed by 
melanization, and eventual metacercariae death (Keeney 
et al. 2007; Koehler and Poulin 2010), but these would 
have to be absorbed by the crabs and disappear from the 
infrapopulation of metacercaria because the percentage 
decreases rather than increases. In that sense, the decline 
in aggregation levels of MZ. nicolli observed in the larg- 
est size-classes of crabs at most locations (Fig. 6A—D) 
would be attributed to the reduction in the abundance of 
metacercariae caused by higher cyst mortality in larger 
crabs. The reductions in mean volumetric metacercariae 
abundance in the largest size-classes at both San Andrés 
and Las Trancas (Fig. 5A, B) resulted from the number 
of large, mature cysts in these crab size-classes without 
overall increase in infrapopulation size. Interestingly, 
host population regulation by parasites is thought be near- 
ly impossible when parasites are highly aggregated in the 
host (Poulin and Morand 2000). 


Conclusions 


Parasite-induced mortality in populations of a marine 
crustacean (E. rathbunae) under field conditions by in- 
frapopulations of the metacercaria of M. nicolli is doc- 
umented in four populations of sand crabs in Guerrero, 
Mexico. These results are consistent with previous stud- 
ies of other Crustacea-Digenea host-parasite systems. 
Evidence that the population of E. rathbunae is affect- 
ed, and possibly regulated, by the infrapopulations of 
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M. nicolli consists of five factors: 1. Mean volumetric 
abundance of metacercariae generally decreased as crabs 
aged (grew larger), but the proportion of mature cysts de- 
creased, suggesting that the metacercariae were growing, 
but that crabs with a higher proportion of more mature 
cysts had been removed (preferentially?) from the pop- 
ulation of crabs; 2. Linear and polynomial regressions 
of mean abundance of helminths (log x+1) vs. ceph- 
alothorax length of crabs were significant for the four 
populations of E£. rathbunae (Table 3); 3. Although both 
the linear and polynomial regressions were highly signif- 
icant, the polynomial regressions had a higher coefficient 
of determination (1) values, indicating that curvilinear 
regressions provided a better fit to the data and confirm- 
ing that selective mortality of hosts was occurring; 4. A 
positive correlation between crab size and abundance of 
helminths was observed at all four locations, confirming 
the potential for host population regulation and/or ma- 
nipulation of host behavior (favoring transmission) by 
parasites; and, 5. The mean number of melanized (dead) 
cysts decreased with the age of the crabs, confirming that 
crabs either were eliminating substantial numbers of par- 
asites throughout their lifetime or more heavily-infected 
crabs were being removed from the populations, but di- 
rect causes of the decrease could not be identified. The 
next steps in this study should be to evaluate the genetic 
composition of the infrapopulations and the metapopula- 
tions of M. nicolli in order to confirm or reject differential 
effects at some locations caused by genetic differences. 
As well, possibility of histopathology should be evaluat- 
ed to determine if it, rather than behavior modification, 
is the cause of the changes in mean abundance between 
size Classes of crabs and the decrease in melanized cysts 
in larger crabs. 
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